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SECTION  I 
INTRODUCTION 

This  manual  is  intended  to  present  the  program  logic  and  user  instruc- 
tions for  a  Data  Storage  Program  (DSP),  with  plotting  and  statistical 
options.   The  Program  was  developed  to  store,  display  and  analyze  both 
raw  and  transformed  data  collected  as  part  of  a  DOT/FHWA  Project, 
Constituents  of  Highway  Runoff  DOT-FH-1 1-8600. 

The  program  allows  researchers  and  other  interested  users  ready  access 
to  the  project  data.   The  program  also  provides  a  readily  available  and 
convenient  means  for  storage  and  analysis  of  environmental  and  operational 
data  collected  during  future  highway  runoff  studies. 

This  manual  is  divided  into  two  sections.   The  first  section  includes  the 
systems  logic,  computer  requirements,  a  description  of  the  variables  and 
the  program  flow  charts.   The  later  section  presents  the  users  manual,   t 
operating  instructions,  and  the  program  source  listing. 

The  manual  has  been  written  to  aid  the  engineer  or  computer  analyst  in 
installing  and  using  DSP  on  their  computer.   The  manual  assumes  that  the 
reader  is  familiar  with  FORTRAN  IV  and  has  a  working  knowledge  of  the  A 
specific  job  control  language  for  the  computer  system  being  utilized. 


SECTION  I  I 
DATA  STORAGE  PROGRAM  (DSP)  SYSTEM 


COMPUTER  REQUIREMENTS 


DSP  is  designed  to  operate  in  a  batch  mode  on  computer  systems  having 
the  following  devices  and  capabilities: 

1 .  A  FORTRAN  compi ler. 

2.  The  ability  to  create  permanent  and/or  temporary  data  sets 
with  random  access  capabilities. 

3-   100K  bytes  of  core  storage  for  DSP. 

k.       180  bytes  of  core  storage  for  DSP  and  plotting  option. 

5-   Access  to  a  peripheral  plotting  device  if  plots  are  desired. 

SYSTEM  OVERVIEW 

DSP  is  designed  to  be  flexible  and.operate  in  various  modes  depending 
on  the  user's  needs.   The  standard  mode  illustrated  in  Figure  1  is  to 
assemble  the  input  on  punched  cards  and  execute  DSP,  storing  the 
plotting  data  and  statistical  data  on  permanent  data  sets.   The  input 
and  output  from  DSP  should  then  be  checked  for  errors.   Once  it  has 
been  determined  that  the  data  is  correct  the  user  may  use  the  plot 
generating  option  illustrated  in  Figure  2.   This  program  option  reads 
the  input  from  the  permanent  data  set  and  generates  the  x,  y  coordinates 
and  pen  movements  for  a  peripheral  plotting  device. 

The  main  program  has  three  subroutines--TIMCON,  FORMUP,  and  F0RMN2 
(Figure  l).   Subroutine  TIMCON  (Figure  3)  is  used  to  convert  military 
time  into  alpha  characters.   Subroutines  FORMUP  and  F0RMN2  are   called 
to  build  the  format  statements  for  each  of  the  thirty  possible  discrete 
or  composite  quality  parameters  that  are  input  to  the  program.   These 
subroutines  therefore  allow  selective  print  format  of  the  desired 
parameters . 

By  using  the  proper  job  control  cards,  permanent  data  sets  may  be 
created  to  store: 

1 .  I nput  card  images  . 

2.  Plotting  information. 

3.  Statistical  information. 

In  this  way,  a  data  handling  and  storage  system  which  can  be  tailored 
to  the  user's  needs  is  built  around  the  DSP.   For  example,  if  the  plot 
generating  option  is  not  utilized  the  JCL  creating  the  plot  files  is 
simply  omitted. 


Call  TlMCON 
Call  FORMUP 
Call  FQRMV 


/JCL  To 
Store  Al I 
Input  As 
Permanent  Data  Set 


«/rite  Line 

Printer 

Quant i  ties 

<Vnd 

IJual  i  ties 

From  Eve 


Write 
i  Plotting  Data 


f 


Write 

Statistical 

Data 
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Figure  1.   DSP  standard  operating  mode. 
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Plot  On 
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Pen  Movements 


Figure   2.      Plot   generating   option    flow   chart 


Subroutine 

TIMCON 
(ITIM,  LIT) 


I 


Dimens  ion 
LIT(i»),  NUM  (10) 


Data  Num 
0-1-2-3-^-5-6-7-8-9 


TIM  o  ITIM 
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TIM  -TIM  *  .0001 


DO  10  I   -1,4 


TIM  =  TIM  *  10 


TIM  ■  TIM  +  .00001 
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ITMP  -  TIM 


TIM  -  TIM  -  ITMP 


DO  10  J  -  1,10 


Yes 


LIT (I)  ■  Num  (J) 


5 


Continue 


A 


f   Return     J 

Figure  3.   Flow  chart  of  subroutine  TIMCON 
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JCL  For  IBM/AMDAHL  370  System 

On  the  370  System  the  Data  Storage  Program  requires  that  data  be 
assigned  to  input  and  output  data  sets  as  follows: 

//FT05F001     is  assigned  to  the  input  data  set,  data  set  5. 
//FT06F001     is  assigned  to  the  output  print  set,  data  set  6. 
//FT09F001     is  assigned  as  an  output  (permanent  or  temporary) 

for  the  plot  generating  program,  data  set  9. 
//FT10F001     is  assigned  to  the  statistical  data  set,  data  set  10, 
//FT40F001     is  assigned  to  the  output  data  set  for  input  to  the 

user's  Plotter,  data  set  ^0. 

The  monitoring  and  background  data  for  each  event,  data  set  5,  is  input 
to  the  data  storage  program.   The  output  print,  data  set  6,  is  generated 
onto  the  line  printer.   The  plotting  program  is  an  optional  step.   If 
the  card  //FT09F001  DD  SYS0UT=A  is  present,  data  set  9  wi 1 1  be  generated 
onto  the  line  printer.   If  instead  the  card  //FT09F001  DD  DSN=name  is 
present,  data  set  9  wi 1 1  be  stored  as  a  permanent  file  (name).   If  the 
card  //FT09F001  DD  DUMMY  or  if  no  card  is  entered  the  plotting  option  is 
cancel  led. 

The  Plot  Generating  program  reads  data  set  9  (see  plot  generating  JCL 
Figure  2)  and  generates  the  actual  x  and  y  coordinates  and  pen  movements. 
The  //FT40F001  DD  DSN=name  JCL  card  directs  the  x  and  y  coordinates  and 
pen  movements  to  be  stored  as  data  set  kO.      The  user's  in-house  plotter 
then  utilizes  data  set  *t0,  FT40F001  DD  DSN=name,  to  create  the  desired 
plots. 

DSP  may  also  be  used  to  generate  statistical  data  for  each  moni tori ng 
event.   The  parameters  available  for  analysis  are  shown  in  Table  1. 
The  statistical  file,  which  is  in  the  form  of  coordinates,  is  written 
on  FT10F001  for  each  event  and  should  be  kept  as  a  permanent  file.   The 
user  may  then  use  the  data  on  the  permanent  files  for  correlation 
analysis  or  to  generate  summary  statistics  (Figure  k) .      Statistical 
programs  such  as  Statistical  Package  for  the  Social  Scientist  (SPSS), 
Statistical  Analysis  System  (SAS)  or  Biomedical  Computer  Programs  (BMD) 
can  be  used  to  perform  the  statistical  analyses. 

DSP,  the  Plotting  Option,  and  the  user's  plotting  software  can  be  linked 
together  as  shown  in  Figure  5tNalthough  this  configuration  is  not  the  ,.\ 
only  possible  scheme  for  setting  up  these  programs. 


Table  1.   Parameters  available  for  statistical  analysis 


Total  precipitation  (cm) 

Average  precipitation  (cm) 

Maximum  intensity  (cm/hr) 

Number  of  dry  days 

Number  of  dry  days  less  than  one  inch  accumulation 

Average  daily  traffic  (vehicles/day) 

Average  traffic  for  period  before  storm  with  less 

than  one  inch  accumulation  (vehicles/day) 

Average  truck  traffic  (trucks/day) 

Air  particulate  fallout  rate  (gm/m2/day) 

Total  runoff  (m3) 

Length  of  event  (minutes) 

Average  runoff  (m3/minute) 

Maximum   runoff    (nwminute) 

Time  into  event  of  maximum  runoff  (minutes) 
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Composite  Hg  (kilograms) 

Composite  SS  (kilograms) 

Composite  VSS  (kilograms) 

Composite  BOD5  (kilograms) 

Composite  T0C  (kilograms) 

Composite  COD  (kilograms) 

Composite  TKN  (kilograms) 

Composite  NO2  (kilograms) 

Composite  NO3  (kilograms) 

Composite  TN  (kilograms) 

Composite  TPO/j  (kilograms) 
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Composite  Zn  (kilograms) 

Composite  Sn  (kilograms) 

Composite  Fe  (kilograms) 

Composite  Cu  (kilograms) 

Composite  Cd  (kilograms) 

Composite  Cr  (kilograms) 
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Composite  PCB's  (kilograms) 
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Figure   /».      Flow  chart   for   using   statistical    programs 
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Figure    5.      Flow  chart   of    DSP   and   PGP    linked    program. 
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JCL  For  Running  and  Storing  DSP  -  The  following  JCL  deck  is  required: 

1.  //JOB  CARD, 

2.  //MESSAGE  CLASS 

3.  //STPNM  EXEC  PGM=D0TPGM,REGI0N=100K,TIME=1 
k.  //FT06F001  SYS0UT=A 

5.  //FT09F001  DD  DSN=PLT01,V0L=SER=NNNNNN,UNIT=DISK, 

SPACE=(CYL,  (1,1)),DCB=(RECFM=FB,LRECL=80, 
BLKSIZE=800),DISP=(NEW,CATLG,DELETE) 

6.  //FT10F001  DD  DSN=STAT1 ,VOL=SER=NNNNNN ,UN IT=D I SK, 

SPACE=(CYL,  (1 ,1)) ,DCB=(RECFM=FB,LRECL=80, 
BLKSIZE=800),DISP=(NEW,CATLG,DELETE) 

7.  //FT05F001  DD  * 

[Tn6eAt  Input  Ccutd  Vnok  fteAe) 

8.  /* 

9.  // 

In  the  above  example,  data  set  9  is  assigned  the  name  PLT01 . 
This  data  set  is  stored  on  disc  which  is  identified  by  NNNNNN  and  has 
an  allocation  of  one  primary  cylinder  and  one  secondary  cylinder.   The 
fixed  block  record  format  has  a  logical  record  length  of  80  and  a  block 
size  of  800.   The  disposition  of  the  data  set  named  PLT01  is: 

1.  It  becomes  a  new  data  set. 

2.  It  will  be  cataloged. 

3.  If  execution  terminates  before  completion  of  the  program, 
the  data  set  will  be  deleted. 

JCL  For  Plot  Generating  Program  -  The  following  JCL  deck  is  required: 

1.  //JOB  CARD, 

2.  //MESSAGE  CLASS 

3.  //STPNM  EXEC  PGM=PCTGEN,REGI0N=100K,TIME=1 

4.  //FT40F001  DD  DSN=PGEN1,V0L=SER=NNNNNN,UNIT=DISC, 

SPACE=(CYL,  (1  ,1)), 

[IvibQAt  Plot  VaIvqjl  VnoQinm  HeAej 

5.  //         DCB(RECFM=FB,LRECL=80,BLKSIZE=800), 

DISP=(NEW,CATLG, DELETE) 

6.  //FT09F001  DD  DSN=PLT01,DISP=(0LD) 
7-   /* 

8.   // 

Data  set  kO    is  stored  as  PGEN1  and  has  the  same  characteristics  as  data 
set  9  from  DSP.   In  the  above  example,  data  set  9  is  an  old  data  set 
being  used  as  input  to  the  plot  generating  program. 


Link  Edit  JCL  For  DSP  -  The  following  JCL  deck  is  required: 

1.  //JOB  CARD, 

2.  //MESSAGE  CLASS 

3.  //STPNM  EXEC  FORTGCLR 
k.  //FORT. SYS  IN  DO  - 

5.  //LKED.SYSIN  DO  * 

6.  ENTRY  MAIN 

7.  NAME  DOTPGM(R) 

This  JCL  deck  is  assembled  to  compile,  link  edit. and  resolve  an  input 
program  using  a  "G"  compiler.   The  entry  point  for  execution  of  the 
program  is  the  "MAIN"  subroutine.   The  program  will  be  stored  as  an 
executable  program  called  DOTPGM  and  will  replace  any  other  program  on 
the  load  library  which  has  the  name  DOTPGM. 

JCL  For  Input  To  The  User's  Plotter  -  The  following  card  is  required: 

//FT40F001  DD  DSN=PGEN1 ,DNP=(0LD) 

JCL  for  Control  Data  7600  System 

Compi le  Step  -  Control  cards  which  make  up  the  compile  portion  are  as 


fol lows: 

1. 

COMP,STTCZ,Pn,Tnnn. 

2. 

ACCOUNT, Tnnnn. 

3- 

REQUEST(LG0,*PF) 

k. 

ATTACH(OLDPL.DSP) 

5. 

UPDATE(F,C) 

6. 

REWIND(COMPILE) 

7. 

FTN(A,I) 

8. 

CATALOG (LG0,DSPLG0) 

9. 

7/8/9 

10. 

6/7/8/9 

Card  one  indicates  the  job's  name,  designation  code,  priority  and  time 
limit.   The  job  in  this  example  (COMP)  is  designated  to  run  on  STTCZ 
(System  Twin  Cities-Z),  which  is  the  identification  code  for  the  CDC 
7600  located  in  Minnesota.   Pn  sets  the  priority  level  where  "n"  can 
have  a  value  of  2  (lowest  priority),  k,    or  6  (highest  priority).   Tnnn 
sets  the  execution  time  (maximum=100  seconds).   Card  two  contains  the 
account  number,  Tnnnn.   The  third  card  reserves  local  file  space  which 
is  designated  LGO.   Upon  completion  of  program  execution  LGO  becomes  a 
permanent  file.   If  execution  terminates  before  program  completion,  the 
data  set  (LGO)  will  be  deleted. 

Cards  four,  five  and  six  must  be  included  if  the  input  data  is  to  be 
obtained  from  a  storage  device.  Card  four  calls  the  data  file  named 
OLDPL  and  the  compiled  program  file  DSP  (Data  Storage  Program)  from  the 

I  I 


storage  device.   Card  five  specifies  that  this  is  a  full  (F)  update, 
rather  than  partial  update.   The  C  indicates  a  compiled  file.   Card 
six  commands  the  device  to  backspace  to  file  start  in  preparation  for 
the  next  program  execution.   Cards  four,  five  and  six  must  be  omitted  if 
cards  are  used  to  input  the  data.   The  input  card  deck  would  be  inserted 
after  card  nine.   The  designation  code  I  on  card  seven  indicates  input 
data.   The  designation  code  A  on  the  same  card  indicates  that  program 
execution  should  be  terminated  if  a  job  control  error  is  detected.   Card 
eight  catalogs  the  data  under  file  name  LGO  and  compile  name  DSPLGO. 
Card  nine  is  an  'end  of  file1  statement,  indicating  that  the  file  is 
complete.   It  is  followed  by  the  program  card  deck  only  if  cards  four, 
five  and  six  were  omitted.   Card  ten  is  an  'end  of  job'  statement,  in- 
dicating that  this  job  is  complete.   The  same  job  can  be  run  for  PLOT, 
changing  DSP  and  DSPLGO  accordingly. 

Execution  Step  -  Control  cards  which  make  up  the  execution  portion  are 


as  fol lows : 

1. 

DSP,STTCZ,Pn,Tnnn. 

2. 

ACCOUNT, Tnnnn. 

3. 

REQUEST(TAPE9,*PF) 

h. 

REQ_UEST(TAPE10,*PF) 

5. 

ATTACH(LG0, DSPLGO) 

6. 

MAP (OFF) 

7. 

LDSET(PRESET=ZERO) 

8. 

LGO. 

9. 

CATALOG (TAPE9,DSP9) 

10. 

CATAL0G(TAPE10,DSP10) 

11. 

7/8/9 

12. 

6/7/8/9 

Cards  three  and  four  create  local  file  spaces  named  TAPE9  and  TAPE10, 
respectively.   After  execution  these  local  files  become  permanent  files 
(*PF).   The  fifth  card  calls  the  data  file  named  LGO  and  the  compiled 
program  file  DSPLGO  from  the  storage  device.   These  files  were  created 
previously  during  the  Compile  Step.   Card  six  terminates  the  MAP  option 
which  would  otherwise  list  relative  storage  locations  of  control  sections 
When  the  map  option  is  used,  an  estimate  is  obtained  of  the  required 
program  region  size.   Card  seven  presets  all  variables  in  core  to  zero. 
Card  eight  executes  LGO.   Card  nine  and  ten  catalog  the  data  under 
file/compile  names  TAPE9/DSP9  and  TAPE10/DSP10,  respectively. 

Plot  Device  Procedure  -  Control  cards  required  for  the  plot  device 
procedure  are  as  follows: 

1.  PLOT,STTCZ,Pn,Tnnn. 

2.  ACCOUNT, Tnnnn. 

3.  REQUEST(TAPEAO,*PF) 

4.  ATTACH (TAPE9,DSP9) 

5.  ATTACH (LGO, PL0TLG0) 
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6.  MAP  (OFF) 

7.  LGO. 

8.  CATALOG (TAPE40,PLOT40) 

Card  three  reserves  a  local  file  space  named  TAPE^O  and,  after  execution, 
it  becomes  a  permanent  file  (*PF) .   Cards  four  and  five  call  the  data 
files/compiled  files  named  TAPE9/DSP9  and  LG0/PL0TLG0,  respectively. 
Card  seven  executes  LGO,  and  card  eight  catalogs  the  data  under  file 
name  TAPE40  and  compile  name  PL0T40. 

Tape  Update  Procedure  -  Control  cards  required  for  tape  update  are: 

1.  TAPE,STTCZ,Pn,Tnnn. 

2.  ACCOUNT, Tnnnn. 

3.  ATTACH (0LDPL1, DSP) 

4.  ATTACH (0LDPL2, PLOT) 

5.  UPDATE(P=0LDPL1,F,N=NEWPL1,W) 

6.  UPDATE(P=0LDPL2,F,N=NEWPL2,W) 

7.  STAGE(TPAE,POST,NT,PE) 

8.  REWIND(NEWPL1,NEWPL2) 

9.  C0PYBF(NEWPL1,TAPE) 

10.  C0PYBF(NEWPL2,TAPE) 

11.  6/7/8/9 

Cards  three  and  four  call  the  data  files/compiled  files  named  0LDPL1/DSP 
and  0LDPL2/PL0T.   Cards  five  and  six  obtain  update  information  from 
data  files  (P)  named  0LDPL1  and  0LDPL2,  respectively.   These  cards  also 
create  new  (N)  updated  local  files  named  NEWPL1  and  NEWPL2,  respectively. 
The  letter  "F"  indicates  that  this  is  a  full  update,  rather  than  a 
partial  update,  while  the  letter  "W"  indicates  that  the  files  are  to  be 
placed  on  the  magnetic  tape  in  sequential  order.   Card  seven  commands 
the  storage  device  to  write  on  a  magnetic  tape  (TAPE),  in  nine  tracks 
(NT),  phase  and  code  at  1 600  bytes  per  inch  (PE) ,  with  even  parity 
(POST).   Card  eight  commands  the  storage  device  to  rewind  the  tape. 
Cards  nine  and  ten  unload  the  data,  which  are  in  local  files  named 
NEWPL1  and  NEWPL2,  onto  a  magnetic  tape  (TAPE)  in  binary  file  language 
(BF). 

Card  Update  Procedure  -  Control  cards  which  make  up  the  update  portion 
are   as  fol lows (cards  to  magnetic  tape): 

1.  TAPE,STTCZ,P6,T100. 

2.  ACCOUNT, Tnnnn. 

3.  UPDATE (F,N=NEWPL1,W) 
k.  UPDATE(F,N=NEWPL2,W) 

5.  STAGE(TAPE,POST,NT,PE) 

6.  REWIND(NEWPL1,NEWPL2) 

7.  C0PYBF(NEWPL1,TAPE) 

8.  C0PYBF(NEWPL2,TAPE) 
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9.      7/8/9 

[In&eAt  Data  Storage.  Vti.oQh.am) 

10.  7/8/9 

[In&ejvt  Vlot  Vn.OQh.am) 

11.  6/7/8/9 

Cards  one  through  eight  were  already  discussed  in  the  previous  section 
on  tape  update  control  cards.   Card  nine  is  an  'end  of  file'  statement. 
In  this  JCL  configuration  it  acts  as  a  partition  between  the  control 
cards  and  the  first  deck  which  updates  DSP.   Card  ten  functions  as  a 
partition  between  the  first  deck  and  the  second  deck  which  update  DSP 
and  PLOT,  respectively. 

DATA  STORAGE  PROGRAM 

Figure  6  is  a  detailed  flow  chart  of  the  data  storage  program  and 
Table  2  lists  the  variable  names  and  their  descriptions.   This  flow 
chart  does  not  include  the  subroutines  or  the  plotting  option.   The 
DSP  is  divided  into  seven  blocks,  with  entry  into  each  block  controlled 
by  a  computed  G0  T0  depending  on  the  value  of  the  variable  KKK. 
Appendix  Table  I  contains  the  Data  Storage  Program  listing  and 
Appendix  Table  I  I  contains  the  plot  program  listing. 

The  first  block,  B1  START  to  B1  END,  is  the  controlling  or  "Executive" 
block.   In  this  block,  all  arrays  are  dimensioned  and  the  first  four 
data  cards  are    read  in.   The  first  card  contains  the  event  number,  a 
snow  melt  option  indicator,  and  a  plot  option  indicator.   The  second 
card  contains  the  area  of  the  site  in  acres,  and  a  20  character  title 
which  appears  on  each  page  of  the  output  and  on  all  plots.   This  title 
normally  contains  site  location,  event  number  and  any  other  pertinent 
information.   The  third  card  contains  the  month,  day  and  year  of  the 
start  of  the  event  and  the  starting  time  of  the  rainfall  and  runoff. 
The  fourth  card  contains  the  date  (month,  day  and  year)  and  time  the 
rainfall  event  ended.   The  remainder  of  the  block  converts  military 
time  into  alpha  characters  using  the  subroutine  TIMCON.   Initially  the 
variable  for  the  computed  G0  T0  statement  (KKK)  is  set  to  a  value  of 
one,  preventing  the  DOT  banner  from  being  written  on  the  first  two 
pages.   The  D0  statement  which  drives  the  program  is  next  (D0  111  I  = 
1,2).   The  remaining  segment  of  program  logic  writes  the  DOT  banner 
(if  KKK  >  1),  run  title,  rainfall  starting  date,  rainfall  starting 
time,  and  event  number  on  each  page  of  output.   This  loop  also 
contains  the  computed  G0  T0  statement  which  determines  which  block  of 
program  logic  is  executed  next. 

Upon  completion  of  the  D0  loop  in  block  one,  control  is  passed  to 
block  two  which  extends  from  B2  START  to  B2  END.   This  block  determines 
if  the  graph  option  will  be  used.   If  the  graph  option  is  selected 
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Figure   6.      Flow  chart  of  data   storage  program. 
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Figure   6.       (continued). 
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Figure   6.       (continued) 
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Figure   6.       (continued) 
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Figure   6.       (continued) 
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Figure   6.       (continued). 
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Table  2.   Variable  names  and  description. 


Variable  name 
ADT 


ADTI 


AREA 
BAREA 


CFSTI 


Fo  rma  t 
Real 


Real 


Descri  ption 


Real 
Real 


CFS  (K) 

Real 

CFST 

Real 

CFSTIM  (K) 

Real 

Real 


CONLB  (JE, 

Real 

IPT,  (JI)) 

CONVRT 

Real 

(JE,  JI) 

CUM 

Real 

DURATN 

Real 

ENDD  (I) 

1  nteger 

ENDT 

1  nteger 

EVENT 


nteger 


Input  variable  for  average  daily  traffic 
since  last  event  (vehicles  per  day) 

Input  variable  for  average  daily  traffic 
for  period  before  the  storm  event  with  an 
accumulation  of  less  than  one  inch  of 
rainfall  (vehicles  per  day) 

Input  variable  for  site  area  (acres) 

One  of  two  possible  values;  1)  area  of 
dustfall  bucket  which  is  read  in,  or 
2)  defaults  to  the  area  of  a  standard 
eight  inch  dustfall  bucket  (0.032^  square 
inches) 

Array  for  storing  input  flows 

Sum  of  the  input  flow  rates  (cubic  feet 
per  second) 

Used  to  associate  the  incremented  flow 
time  with  the  appropriate  input  flow 

Total  flow  (cubic  feet)  as  defined  by  the 
equation:  CFST  *  time  interval  *  60 

Array  for  storing  discrete  quality 
pollutant  loading  values  (lbs/min) 

Array  for  computing  discrete  quality 
pollutant  loading  values  (lbs/min) 

Cumulative  rainfall  sum  (inches) 

Input  variable  for  duration  of  the 
rainfall  event  (hundredths  of  an  hour) 

Three  position  array  for  storing  end  of 
event  (month,  day,  year) 

Input  variable  for  time  event  ended 
(military  time,  e.g.,  2^00  hours) 

Input  variable  for  event  number 
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Table  2.   (continued) 


Variable  name 
FLOW  (I) 


FLOWST 


IADT1 


Format 
Real 


Rea 


Descri  pt ion 


GALLON 

Real 

GMSAVE 

Real 

GMSSQ 

Real 

GMS1 

Real 

GMS2 

Real 

GMS3 

Real 

HYETO 

Alpha- 

Numeri  c 

I 

1  nteger 

IADT 

Integer 

I  nteger 


ICOL 

Integer 

ICOUNT 

1  nteger 

IEOD 

1 nteger 

IEVENT 

Integer 

IHOUR 

1  nteger 

IIII 

Integer 

Array  for  storing  input  flow  values  (cubic 
feet  per  second) 

Input  variable  for  flow  start  time 
(military  time)  from  which  all  other  flow 
times  are  calculated  using  the  input  flow 
time  interval 

Total  flow  in  gallons 

Average  grams  of  dust  for  number  of 
buckets  sampled 

Air  particulate  fallout  rate  (gm/m^/day) 

Grams  of  dust  for  bucket  #1 

Grams  of  dust  for  bucket  #2 

Grams  of  dust  for  bucket  #3 

If  IOPT  =  HYETO,  the  hyetograph  data  will 
be  written  to  the  plotting  file 

Array  positional  value  (row) 

Rounded  up  or  down  integer  value  of  ADT 
which  is  read  in  as  a  real  number 

Rounded  up  or  down  integer  value  of  ADT1 
which  is  read  in  as  a  real  number 

Array  positional  value 

Counts  number  of  input  rainfall  values, 
used  to  compute  average  rainfall 

Used  as  end  of  data  set  (-999) 

Integer  value  of  real  variable  "EVENT" 
used  for  plotting 

Time  in  hundreds  of  hours,  e.g.,  2200  =  22 

Array  positional  value 
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Table  2.   (continued) 


Variable  name 


Format 


IINT 

Integer 

IMINUT 

1 nteger 

INPUTQ  (JE) 

Alpha- 

Numeric 

INQSYM 

Alpha- 

Numeric 

Description 


IOPT 


ISTHR 


LB  (JE) 

LNUM, 

LNUM1,  and 
LNUM2 


Integer 


nteger 


ISTMIN 

I  nteger 

ISWTCH 

Integer 

JE 

I  nteger 

K 

I  nteger 

KG  (JE) 

Integer 

KI 

Integer 

KJ 

Integer 

KKK 

Integer 

Real 

Alpha- 
Numer  i  c 


Integer  value  of  RINTVL,  rainfall  interval 

Used  to  save  integer  value  for  minutes 
part  of  military  time  designation,  nnnn 

Input  variable  for  discrete  quality 
abbreviations,  e.g.,  SS  FE,  COD,  etc. 

Data  statement  containing  30  possible 
input  quality  parameters  and  k   blanks 

Graph  option  if  IOPT  =  IGRAPH  the  program 
writes  plotting  data  on  Unit  9  for  IHYETO, 
IHYDRO,  and  IPLUTO 

Input  variable  for  hour  rainfall  started 
(mi  1 i  tary  time) 

Input  variable  for  minute  rainfall  started 

Position  pointer,  counter 

Array  positional  value 

Value  set  equal  to  zero  and  incremented 
by  one  for  determining  the  first  time 
through  a  loop;  later  used  as  the  length 
of  another  DO  loop. 

Composite  quality  (kilograms) 
Positioner  for  various  arrays 

Positioner  for  various  arrays 

KKK  is  used  to  transfer  from  one  section 
of  program  to  another  depending  upon  its 
val  ue 

Composite  quality  (pounds) 

Transfer  variables  for  the  subroutine 
TIMCON  representing  the  alpha  equivalence 
of  time;  these  variables  are  also  used 
for  plotting 
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Table  2.   (continued) 


Variable  name 
LOC 


Format 

Alpha- 
N  ume  r  i  c 


Descri  pt ion 


Input  variable  for  two  letter  optional 
site  abbreviation 


MN,  M1 
M2  &  M3 


Integer 


M3 

Real 

M3T 

Real 

MGLD1 

Real 

MGLD2 


MGLD3 


MILSC1 


MILSC2 


MILSC3 


Real 


Real 


Real 


Real 


Real 


Main  body  program  pointers  and  counters 
used  to  hold  program  position  while  the 
subroutines  FORMUP  and  F0RMN2  which  build 
formats  for  the  quality  parameters  are 
cal led 

Flow  in  cubic  feet  per  minute 

Total  flow  in  cubic  meters  (M3T  *  IINT) 

Input  variable  for  milligrams  per  liter  of 
dust  and  dirt  in  bucket  No.  1 

Input  variable  for  milligrams  per  liter  of 
dust  and  dirt  in  bucket  No.  2 

Input  variable  for  milligrams  per  liter  of 
dust  and  dirt  in  bucket  No.  3 

Input  variable  for  milliliters  of  samples 
collected  from  dustfall  bucket  No.  1 

Input  variable  for  milliliters  of  samples 
collected  from  dustfall  bucket  No.  2 

Input  variable  for  milliliters  of  samples 
collected  from  dustfall  bucket  No.  3 


NB 


NUMBD 


NUMDD 


NUMDL1 


Real       Input  variable  for  number  of  dustfall 

bucket  samples  associated  with  this  event 

Real       Input  variable  for  number  of  days  the 

dustfall  buckets  were  up  collecting  the 
sample 

Integer     Input  variable  for  number  of  dry  days 
before  this  event 

Integer     Input  variable  for  number  of  dry  days 

before  this  event  with  less  than  one  inch 
rainfall  accumulated 

Real       Cumulative  rainfall  sum  (centimeters) 
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Table  2.   (continued). 


Variable  name 


Format 


PLUTO 

Alpha- 

Numeric 

QUALTY(I,M) 

Real 

RAIN (I) 

Real 

RAIN1  (I) 

Real 

RAIN1  (II) 

Real 

Description 


RINTVL 

RMAX 
RSUM 


RTIM1 
RTIM2 

RTIM3 

RTIME 

RTITLE 

SAMPLE (I) 

SNOW 


Real 

Real 
Real 


Real 
Real 

Real 

Real 

Alpha- 
Numeric 

Integer 
I nteger 


If  IOPT  =  PLUTO,  pol lutographs  will  be 
plotted 

Array  for  storage  of  quality  parameters 

Array  for  storing  rainfall  record 

Used  to  hold  the  absolute  value  of  the 
differences  between  the  rainfalls 

Rainfall  value  for  time  one  minus  rainfall 
value  for  time  two 

Input  rainfall  interval  for  rainfall 
record  (minutes) 

Maximum  rainfall  intensity  for  an  event 

Variable  used  in  two  different  program 
sect  ions : 

1)  Rainfall  Sum  (centimeters) 

2)  Flow  Sum  (cubic  meters) 

Actual  flow  start  time 

Time  of  maximum  flow  (hours,  minutes); 
later  converted  to  minutes 

Time  of  last  flow  in  the  loop  (hours, 
minutes) 

Current  time  used  to  convert  flow  start 
time  in  the  cases  where  IHOUR  >  2k   or 
IMINUT  >  60 

Site  name 


Input  variable  for  sample  number  (usually 
a  laboratory  identification  number) 

Input  variable  with  two  possible  values: 

1)  One  signals  data  from  snow  melt 

2)  Blank  signals  rainfall  event 
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Table  2.   (continued) 


Variable  name 
STARTD(I) 


Format 
I nteger 


Description 


STARTD 

I  nteger 

STARTT 

Integer 

SUM 

Real 

T1 

Real 

T2 

Real 

T3 

Real 

TEMP  (I) 

Real 

TIME  (I) 

Integer 

TINT 

I  nteger 

TINT2 


TMP 


TMP1 


TMP2 


Integer 


Real 


Real 


Rea 


TOTRN 

Real 

TPRCP 

Real 

TPRECP 

Real 

TRUCKS 

Real 

XM3 


Real 


Three  position  array  for  event  start 
(month,  day,  year) 

Input  variable  for  overflow  start  time 

Input  variable  for  rainfall  start  time 

Variable  used  to  sum  rainfall 

RTIM1  in  minutes 

RTIM2  in  minutes 

RTIM3  in  minutes 

Array  for  storing  rainfall  intensity 
values  (inches  per  hour) 

Array  to  store  time  values 

Time  interval  between  discrete  samples 
for  sampler  No.  1 

Time  interval  between  discrete  samples 
for  sampler  No.  2 

Storage  location  for  first  rainfall  input 
value  used  for  determining  length  of  event 

The  running  value  of  rainfall  intensity 

in  centimeters;  used  in  a  loop  to  determine 

maximum  rainfall  intensity 

The  running  value  of  the  absolute  value  of 
the  difference  between  rainfall  value  in 
cent  imeters 

Input  variable  for  total  rainfall  (inches) 

Total  rainfall  (centimeters) 

Total  rainfall  volume  (cubic  feet) 

Input  variable  for  average  truck  traffic 
(trucks  per  day) 

Flow  (cubic  feet  per  minute) 
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(IOPT  =  IGRAPH) ,  the  program  will  write  to  the  line  printer  that  the 
plotting  option  is  in  effect.   This  block  of  program  logic  will  also 
write  the  title  of  the  run,  event  number,  event  starting  date  and 
rainfall  start  time  to  a  permanent  data  file  which  has  been  set  up  by 
the  user's  JCL  (Figure  2).   If  the  graph  option  is  not  selected 
(IOPT  ±   IGRAPH),  control  will  bypass  the  plot  write  statements.   KKK  is 
set  equal  to  2  and  control  is  returned  to  statement  107. 

The  computed  G0  T0  statement  then  moves  execution  to  B3  START.   The 
program  logic  of  block  three  writes  the  event  starting  date,  event 
ending  date,  rainfall  start  time,  rainfall  end  time  and  flow  start  time 
to  the  line  printer.   Card  number  five  is  also  read  and  contains  the 
following  information: 

1.  Total  rainfall  volume  (inches). 

2.  Rainfall  duration  (hundredths  of  an  hour). 

3.  Number  of  dry  days  before  this  event  (days). 

^t.   Number  of  dry  days  before  this  event  with  an  accumulative 
rainfall  less  than  one  inch  (days). 

5.  The  time  interval  between  flow  samples  for  sampler  number 
one  (minutes) . 

6.  The  time  interval  between  flow  samples  for  sample  number 
two  (minutes) . 

7-   The  Average  Daily  Traffic  (vehicles  per  day). 

8.  The  Average  Daily  Traffic  during  the  dry  period  with  less 
than  one  inch  of  accumulated  rainfall  (vehicles  per  day). 

9.  Average  truck  traffic  (trucks  per  day). 

The  next  two  IF  statements  are  used  to  round  up  to  a  whole  number  the 
average  daily  traffic  and  average  daily  traffic  during  the  dry  period 
with  less  than  one  inch  of  accumulated  rainfall.   Next,  total  rainfall 
in  inches  is  converted  to  total  rainfall  in  centimeters  and  the  total 
rainfall  volume  in  cubic  feet  is  computed.   After  these  calculations 
are  completed,  the  drainage  area,  appropriate  precipitation  and  dry 
day  data,  and  flow  sampler  time  intervals  are  written  to  the  line 
printer.   Next,  precipitation  data,  dry  day  data  and  traffic  data  are 
stored  in  the  statistical  array  (STAT).   KKK  is  set  equal  to  3  and 
control  is  returned  to  statement  107- 

The  next  block  (Bk    START  to  B^t  END)  contains  the  equations  and  conver- 
sion factors  required  to  calculate  the  air  particulate  fallout  rate. 
First  the  dustfall  banner  is  written  to  the  line  printer.   Next,  card 
number  six  is  read  which  contains  number  of  dustfall  buckets, 
number  of  days  these  buckets  were  in  place,  and  area  of  dustfall  bucket 
in  square  meters,  if  different  from  the  standard  eight  inch  dustfall 
bucket.   A  test  is  made  to  determine  the  number  of  dustfall  buckets 
that  are  to  be  read  in.   If  the  number  is  less  than  or  equal  to  zero 
the  series  of  dustfall  calculations  are  skipped.   If  the  number  of 
dustfall  buckets  is  greater  than  zero,  dustfall  data  for  each  dustfall 
bucket  will  then  be  read  in,  one  card  for  each  bucket  containing  the 
milliliters  of  sample  collected  and  milligrams  per  liter  of  dust  and 
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dirt  for  the  associated  sample.   The  program  then  converts  this  data 
into  grams  of  dust  and  dirt  as  follows: 

Grams  of  dust  and  dirt  =  sample  volume  (ml)  '-'■'  concentration  (mg/1) 

1,000,000 

The  grams  of  dust  and  dirt  collected  by  each  of  the  dustfall  buckets 
at  a  particular  site  is  then  averaged: 

Average  weight  (grams)  = 

sample  weight  (bucket  no.  1)  +  sample  weight  (bucket  no.  2)  + 

sample  weight  (bucket  no.  3) 

number  of  dustfall  buckets  analyzed 

Program  logic  assigns  a  zero  weight  to  buckets  for  which  no  dustfall 
data  was  entered  by  setting  all  MGLD  variables  (mg/1  dust  and  dirt)  to 
zero  before  the  calculations  are    initiated.   The  program  then  writes 
to  the  line  printer  the  appropriate  dustfall  data  which  is  read  in 
for  each  dustfall  bucket  and  the  grams  of  dust  and  dirt  associated 
with  each  sample  along  with  the  number  of  days  the  dustfall  buckets 
were  in  place  collecting  data.   The  air  particulate  fallout  rate  is 
then  computed  for  the  collection  period; 

Air  particulate  fallout  rate  (gm/m^/day)  = 

Average  weight  (grams)/bucket  area    (rrr) 

Number  of  days  dustfall  buckets  were  in  place 

Note:   The  value  0.0324  is  the  surface  area   of  a  standard  8  inch 
dustfall  bucket  in  square  meters  (default  value). 

The  program  then  writes  to  the  line  printer  the  air  particulate  fallout 
rate,  and  stores  the  value  in  the  statistical  array  (STAT).  KKK  is  set 
equal  to  k   and  control  is  transferred  to  statement  number  107. 

The  next  block  of  program  logic  controls  the  rainfall  analysis  (B5 
START  to  B5  END).   This  block  first  reads  card  number  eleven  which 
contains  the  rainfall  interval  (time  in  minutes  between  rainfall  input 
values),  rainfall  start  time  and  hyetograph  plotting  option  (IHYET0) . 
Upon  completion  of  this  read  statement,  the  rainfall  banner  and 
rainfall  interval  are  written  to  the  line  printer.   The  next  step  is 
to  read  the  rainfall  data  (inches)  from  the  rainfall  cards.   Each 
card  contains  eight  rainfall  values  RAIN  (i),  I  =  1,8.   The  last 
value  in  the  rainfall  record  must  end  with  a  -99  to  signal  the  end 
of  the  read  statement. 

A  time  value  is  associated  with  each  rainfall  value  read  in,  using 
rainfall  start  time  and  rainfall  interval.   The  running  sum  (SUM)  and 
accumulated  sum  CUM  (l)  are  associated  in  D0  loop  200.   The  rainfall 
record  is  then  converted  to  metric  (inches  *  2.5*0  and  total 

29 


rainfall  and  maximum  rainfall  intensity  are  found.   The  rainfall  record 
with  associated  times  is  then  written  to  the  line  printer.   If  the 
hyetograph  option  is  called  for  (IOPT  =  IHYETO) ,  the  information  is 
also  written  to  a  storage  file.   The  maximum  rainfall  and  average 
intensity  are   written  to  the  statistical  array  and  control  returns  to 
statement  107- 

The  next  block  of  program  logic  deals  with  the  flow  data  (B6  START  to 
B6  END).   First  the  flow  banner  is  written  out  to  the  line  printer. 
Next,  card  number  thirteen  which  contains  the  flow  start  time  (hours),  flow 
interval  (minutes  between  values),  and  hydrograph  plotting  option  is 
read.   If  IOPT  =  IHYDR0  the  input  hydrograph  information  will  be  written 
to  a  storage  file.   Next,  the  flow  record  is  read  in  one  value  at  a  time 
with  a  maximum  of  eight  values  per  card.   The  program  will  continue  to 
read  these  cards  until  it  encounters  a  ~99-  on  the  last  flow  card. 
When  this  is  encountered,  control  shifts  to  writing  flow  summaries. 
Actual  flow  values,  FL0W  ( I)  are    read  in  a  list  directed  read 
statement  with  the  array  positioner  "I"  ranging  from  one  to  eight  for 
each  value  on  a  specific  card.   These  flow  values  are  then  moved  to 
another  array  CFS  (K) ;  here  the  K  array  positioner  is  assigned  a 
value  from  one  to  the  number  of  flow  points  which  are    read    into  the 
program.   Array  positioner  K  is  incremented  each  time  a  flow  value 
is  processed.   Each  flow  value  CFS  (K)  is  converted  from  cubic  feet 
per  second  (cfs)  to  cubic  meters  per  minute  by  multiplying  CFS  *  1.699- 
The  cubic  feet  per  minute  values  are  then  summed.   The  cfs  values  are 
converted  to  total  gallons  as  follows: 

Total  gallons  =  Zcfs  «  time  i  nterva  1  (sec)  *  4^8.8/7-^8  gal   ...  60  sec 

\   f t3         min 

Each  flow  value  is  assigned  a  time  by  using  the  flow  start  time  and 
incrementing  it  by  the  flow  interval  for  each  value. 

After  total  cubic  feet  (CFSTl) ,  total  cubic  meters  (M3T) ,  and  total 
gallons  (GALL0N)  have  been  computed,  a  flow  summary  of  these  values  is 
written  to  the  line  printer.   Next,  a  maximum  flow  in  cubic  meters  and 
an  average  flow  in  cubic  meters  is  calculated  and  stored  in  the  statis- 
tical array  (STAT).   The  time  variable  flow  data  is  then  written  to  the 
line  printer  in  cubic  feet  per  second  and  cubic  meters  per  minute.   The 
values  for  total  cubic  meters  (M3T) ,  length  of  the  event  in  minutes, 
average  flow  in  cubic  meters,  maximum  flow  in  cubic  meters  (RMAX)  and 
time  in  minutes  from  start  of  storm  to  maximum  flow  are   written  to  the 
statistical  array  (STAT). 

Next,  the  discrete  flow  quality  cards  are  read  in.   A  blank  card  must 
be  inserted  into  the  run  deck  if  no  discrete  quality  parameters  are 
to  be  read  into  the  program.   The  blank  card  causes  program  logic 
to  increment  KKK  to  7  and  transfer  control  back  to  statement  107  to 
read  in  composite  quality.   If  discrete  quality  data  is  available, 
the  program  reads  five  discrete  quality  sample  cards  per  sample  time. 
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The  first  card  contains  the  discrete  analysis  time,  laboratory  sample 
number  and  the  first  seven  of  thirty  possible  discrete  quality  parameters 
The  second  and  third  cards  contain  eight  possible  discrete  flow 
quality  parameters  while  seven  possible  parameters  can  be  listed  on 
the  fourth  card.   This  procedure  is  repeated  for  each  discrete  sample 
analysis  time  until  the  program  reads  a  -99.  at  which  point  the  program 
will  G0  T0  100.   Discrete  flow  quality  values  are   stored  in  the  array 
QUALTY  (I,M)  where  I  indicates  the  discrete  sample  time  and  the  M's 
are    the  appropriate  sample  values. 

It  is  possible  that  some  discrete  analyses  were  performed  for  one  time 
period  and  not  another.   When  program  logic  detects  this  condition,  it 
sets  that  parameter  to  the  value  of  1.0  *  10°0.   Because  this  value  is 
much  too  large  for  the  format  statement  defining  discrete  quality  para- 
meter values,  astericks  appear  in  that  position  for  that  time  period. 
A  program  note  appears  under  the  quality  heading  on  the  output  to 
explain  this  situation. 

The  program  then  calls  the  subroutine  FORMUP  which  builds  the  format 
statements  needed  for  the  selected  discrete  flow  quality  parameters 
used.   The  discrete  quality  parameters  analyzed  are  then  converted 
to  pounds  per  minute  as  follows: 

lb/min  =  discrete  quality  concentrations  (mg/1)  * 

the  appropriate  flow  rate  (cfs)  at  that  time  *  0.0037^6 

(0. 0037^6  converts  mg/1  *  cfs  to  lb/min) 

For  each  discrete  quality  parameter,  the  program  writes  out  sample  time, 
identification  number,  parameter  analyzed  and  its  value  in  both  mg/1 
and  lb/min. 

The  program  next  reads  the  pollutograph  option,  card  number  20.   If 
the  option  I0PT  equals  IPLUT0  the  program  writes  the  pollutograph 
data  to  storage,  and  reads  cards  21  and  22.   These  two  cards  are  used 
to  designate  the  abbreviations  for  the  selected  discrete  quality 
parameters  which  will  appear  on  the  pol 1 utographs .   The  program  then 
writes  mg/1  and  lb/min  data  to  be  plotted  for  all  parameters  analyzed. 
The  variable  KKK  is  set  equal  to  7  and  control  is  transferred  back 
to  statement  number  107- 

The  next  block  (B7  START  to  B7  END)  handles  the  flow  composite  quality 
data.   The  program  first  writes  out  the  composite  quality  banner  and 
then  reads  in  card  number  23  containing  the  composite  laboratory  sample 
number  and  the  first  seven  composite  quality  parameters.  There  must  be 
a  total  of  five  cards  as  the  program  will  continue  to  read  composite 
quality  until  it  reads  a  -99 •  on  the  fifth  card.  The  program  then 
writes  the  composite  laboratory  sample  number  and  reads  the  remaining 
four  composite  quality  cards  (similar  to  discrete  quality  cards).   The 
second  card  contains  the  appropriate  8th  through  15th  sample  values, 
the  third  the  1 6th  through  23rd,  the  fourth  the  2^th  through  the  30th, 
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and  the  fifth  the  -99-  value.   The  program  then  begins  building  the 
format  statements  needed  to  write  out  the  appropriate  flow  composite 
quality  parameters.   The  composite  parameters  are  then  converted  to 
kilograms  (kg)  as  follows: 

Sample  value  in  kg  =  total  cubic  feet  discharged  * 

28.32  (liters  per  cubic  feet)  *   sample  concentration  in  mg/1 

(1  *  10"6) 

and  then  to  pounds  (lbs)  as  follows: 

_    ,    i-i.    sample  value  in  kg 

Sample  value  in  lb  =  £  ,  r_,  r — ; — 

K  0.453^5  /kgv 

Mb' 

The  appropriate  composite  samples  are  then  written  to  the  line  printer 
in  mg/1,  pounds  and  kilograms.   The  program  then  stores  all  composite 
quality  parameters  in  kilograms  in  the  statistical  array  (STAT)  and 
checks  the  option  (IPLUT0) .   If  the  option  is  not  equal  to  IPLUT0  the 
program  writes  a  -999  to  the  storage  file  for  pollutograph  plotting. 
The  -999  indicates  to  the  plot  generating  program  that  this  is  the 
end  of  the  plot  data.   The  program  then  writes  to  the  line  printer 
END  OF  DOT  and  program  flow  for  DSP  is  halted. 
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SECTION  I  I  I 

DATA  STORAGE  PROGRAM  (DSP) 
USER'S  MANUAL 


The  Data  Storage  Program  was  developed  as  part  of  a  FHWA/DOT  Research 
and  Development  project  to  determine  the  characteristics  of  highway 
runoff.   Due  to  the  nature  of  the  study,  an  extensive  data  collection 
and  chemical  analysis  program  was  required  and  more  than  1^,000  pieces 
of  data  were  generated.   DSP  was  developed  to: 

a.  Retain  a  large  quantity  of  field  generated  data  in  an 
organized  and  easily  accessible  format. 

b.  Permit  conversion  of  standard  units  to  metric. 

c.  Provide  a  mechanism  by  which  the  data  can  be  recalled  for 
optional  plotting  or  statistical  analysis. 

The  data  base  developed  during  the  project  is  perhaps  the  most  extensive 
highway  runoff  quantity/quality  data  in  existence.   This  data  is 
available  to  interested  parties  on  magnetic  tape  with  access  provided 
by  DSP.   It  is  further  thought  that  DSP  will  be  valuable  as  a  vehicle 
for  data  storage  and  retrival  to  anyone  embarking  on  highway  runoff 
studies. 

In  its  basic  form  DSP  is  merely  a  convenient  method  for  collating  and 
producing  a  report  of  data  collected  in  a  highway  runoff  study.   Storage 
of  data  by  use  of  DSP  will  permit  a  user  to  apply  optional  plotting 
and/or  statistical  techniques  on  individual  runoff  events  or  series  of 
events . 

Highway  runoff  is  highly  variable  in  that  it  results  primarily  from  pre- 
cipitation or  snowmelts  and  as  such  its  description  contains  a  rainfall 
hyetograph,  a  runoff  hydrograph,  various  quality  parameters  (described 
by  either  discrete  or  composite  quality  analyses),  dustfall  data  and 
various  site  characteristics.   A  description  of  how  to  set  up  input  for 
DSP  follows. 

DESCRIPTION  OF  THE  INPUT 

Input  to  the  data  storage  program  is  coded  onto  27  card  groups,  each 
consisting  of  one  or  more  cards.   Table  3  is  an  abbreviated  listing  of 
these  card  groups  for  quick  reference.   Table  h   contains  a  more  detailed 
explanation  including  format  requirements,  card  columns  used  and 
corresponding  variable  names. 

Several  of  the  values  are  inputted  more  than  once.   For  example,  "start 
of  rainfall"  is  coded  onto  card  group  three  and  card  group  eleven. 
Table  5  lists  a  few  guidelines  for  setting  up  the  data  storage  system. 
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Table  3-   Summary  of  card  groups. 


Card  group  Description  of  card  groups 

1.  Event  number,  smowmelt  indicator  and  plot  option 

2.  Drainage  area  and  title 

3-  Starting  date,  start  of  rainfall  and  start  of  overflow 

4.  Ending  date  and  end  of  rainfall 

5.  Blank  card  for  additional  parameters 

6.  Total  rainfall,  rain  duration,  dry  days,  discrete  time 
interval  and  average  daily  traffic 

7.  Number  of  dustfall  buckets,  days  of  sample  and  area  of 
bucket 

8.  Milliliters  of  sample  and  milligrams  per  liter  in 
bucket  1 

9.  Milliliters  of  sample  and  milligrams  per  liter  in 
bucket  2 

10.  Milliliters  of  sample  and  milligrams  per  liter  in 
bucket  3 

11.  Rainfall  interval,  start  of  rainfall  and  hyetograph 
option 

12.  Rainfall  data  at  each  interval 

13-         Start  of  overflow,  flow  interval,  and  hydrograph  option 

14.         Flow  data  at  each  interval 

15-         Discrete  sample  time,  laboratory  number,  coliform,  pH, 
HG,  SS,  VSS,  B0D5,  TOC 

16.  COD,  TKN,  NO  ,  NO   TN ,  TPO  ,  CL",  PB 

17-  ZN,  SN,  FE,  CU,  CD,  CR,  TS ,  0&G 

18.  PCB,  Pesticides  and  flow  data  for  discrete 

19-  Terminates  discrete  input 

20.  Pollutograph  option 

21.  Additional  plot  information 

22.  Additional  plot  information 

23.  Composite  sample  number,  coliform,  pH ,  HG,  SS,  VSS, 
BOD   TOC 

2k.  COD,  TKN,  N02,  NO   TN ,  TPO  ,  Cl_",  PB 

25.  ZN,  SN,  FE,  CU,  CD,  CR,  TS ,  0&G 

26.  PCB  and  Pesticides 

27.  Terminates  composite  input  and  signals  end  of  job 
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Table  4.   DSP  input  format. 


Card 

Card 

group 

Fo  rma  t 

col umns 

1 

A2 

1-2 

14 

5-8 

12 

14-15 

Descri  pt ion 


A5 


14 

14 
312 

14 

F10.2 
F10.2 

F10.2 
F10.2 

15 


22-26 


10.0 

1-10 

5A4 

21-40 

312 

1-6 

9-12 

15-18 
1-6 

9-12 

1-10 
11-20 

21-30 
31-40 


Site  location,  two  letter 
abbreviation  of  site 

Storm  event  number 

Snowmelt  indicator  (-1  if  event  was 
a  s  nowme 1 1 ) 

Graph  option  indicator  ('GRAPH'  if 
option  is  desired,  otherwise  leave 
blank) 

Area  of  site  (acres) 

Title  i.e.  Los  Angeles  HY-5 

Date  storm  began  (month,  day,  year) 
i.e.  091576 

Rainfall  start  time  (military) 
i.e.  0800 

Time  overflow  started 

Date  storm  stopped  (month,  day, 
year)  i.e.  091676 

Time  storm  stopped  (military  time) 
i .e.  0200 

Total  rainfall  (inches) 

Duration  of  storm  (to  the  nearest 
1/100  hour)  i.e.  1.25  =  1  hour 
15  minutes 

Number  of  dry  days  before  the  storm 

Number  of  dry  days  before  the  storm 
with  an  accumulation  of  less  than 
one  inch  of  rainfall 


41-45   Not  used 


Variable 
name 

L0C 


EVENT 
SNOW 

I0PT 

AREA 
RTITLE 
STARTD 

STARTT 

START0 
ENDD 

ENDT 

T0TRN 
DURATN 

NUMDD 
NUMDI 

NUMDS 
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Table  4.   (continued). 


Card 
group 


Format 

15 

15 

15 
no. 2 

no. 2 
no. 2 
no.o 

F10.0 
F10.0 

F10.2 


F10.2 


9     F10.2 


Card 
col umns 

46-50 


51-55 

56-60 
61-70 

71-80 
1-10 
1-10 

11-20 
21-30 

1-10 


Descr i  pt  ion 


11-20 


1-10 


F10.2 


11-20 


Time  interval  between  discrete 
samples  (in  minutes)  for  sampler 
number  one 

Time  interval  between  discrete 
samples  (in  minutes)  for  sampler 
number  two 

Not  used 

Average  traffic/dry  day  since 
last  storm  event 

Average  traffic  for  period  before 
the  storm  with  an  accumulation  of 
less  than  one  inch  of  rainfall 

Insert  a  blank  card  (must  be 
included).   Reserved  for  additional 
data  or  parameters 

This  card  must  be  included  even  if 
blank.  Number  of  dustfall  buckets 
to  be  analyzed  (Maximum  =  3) 

Number  of  days  buckets  were  up 

Area  of  dustfall  bucket  if  not 
standard 

(if  NB  >  0  (card  group  6),  include 
card  group  7)-   Milliliters  of 
sample  collected  for  dustfall 
bucket  number  one 

Milligrams/liter  of  dust  from 
dustfall  bucket  number  one 

(if  NB  >  1  (card  group  6),  include 
card  group  8).   Milliliters  of 
sample  collected  for  dustfall 
bucket  number  two 

Milligrams/liter  of  dust  from 
dustfall  bucket  number  two 


Variable 

name 

TINT 


TINT2 

NUMPMS 
ADT 

AST  I 


NB 

NUMBD 
BAREA 

MILSCI 


MGLDI 


MILSC2 


MGLD2 
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Table  k.       (continued) 


Card 

group   Format 

10     F10.2 


11 


12 


13 


F10. 2 

F10. 2 

212 

A5 

F10.2 


F10.2 


F10.2 


F10.2 


F10.2 


F10.2 


Card 
col umns 

1-10 


Descr ipt  ion 


11-20 


1-10 


22-26 


1-10 


11-20 
21-30 
71-80 


Ik 

\-k 

Ik 

7-10 

A5 

22-26 

(if  NB  >  2  (card  group  6),  include 
card  group  9)-   Milliliters  of 
sample  collected  for  dustfall 
bucket  number  three 

Milligrams/liter  of  dust  from 
dustfall  bucket  number  three 

Rainfall  interval  (in  minutes) 
for  the  rainfall  record 


17-20    Time  of  rainfall  start,  i.e.  0800 


Hyetograph  plot  option  indicator 
('HYETO'  if  plot  is  desired) 

Rainfall  measurement  for  the  first 
time  interval  (actual  inches  from 
the  rainfall  chart,  this  is  NOT 
inches/hour;  inches/hour  will  be 
computed) 

Rainfall  measurement  for  the 
second  time  interval 

Rainfall  measurement  for  the 
thi  rd  t  ime  interval 

Rainfall  measurement  for  the 
eighth  time  interval 

Rainfall  measurement  for  the 
Nth  time  interval 

-99- ,  this  value  terminates  the 
rainfall  record,  it  may  occur  in 
any  field  and  on  any  card 

Flow  start  time  (military  time) 

Time  interval  between  flow  values 

Hydrograph  plot  option  indicator 
(■HYDRO1  if  plot  is  desired) 


Variable 
name 

MILSC3 


MGLD3 


RINTVL 


ISTHR  & 
ISTMIN 

I  OPT 


RA I  N  ( 1 ) 


RAIN(2) 
RAIN(3) 
RAIN(8) 
RAIN(N) 
RA!N(N+1) 

FL0WST 
INT 
I  OPT 


37 


Table  4.   (continued) 


Card           Card  Variable 

group  Format   col  umns   Descr  ipt  ion name 

14  F10.2      1-10    Flow  value  (in  CFS)  for  the  first     FL0W(1) 

time  interval 

F10.2     11-20    Flow  value  (in  CFS)  for  the  second    FL0W(2) 

time  interval 

F10.2     21-30    Flow  value  (in  CFS)  for  the  third     FL0W(3) 

time  interval 

F10.2     71-80    Flow  value  (in  CFS)  for  the  eighth    FL0W(8) 

time  interval 

F10.2  Flow  value  (in  CFS)  for  the  Nth       FLOW(N) 

time  interval 

F10.2  -99-,  This  value  terminates  the      FL0W(N+1) 

flow  record,  it  may  occur  in  any 
field  and  on  any  card 

15  14       1-4     Time  first  discrete  quality  sample    T I  ME ( 1 ) 

was  taken  (military  time)  i.e.  0915 

l4       7-10    Sample  ID  no.  for  discrete  quality    SAMPLE(1) 
samples  (usual ly  an  ID  number 
assigned  by  the  laboratory  performing 
the  quality  analyses) 

E10.4     11-20    Coliform  value  for  the  first        QUALTY(1,1) 

sample  (C0LF)  (number/100  ml ) 

F10.2     21-30    pH  value  for  the  first  sample(PH)    QUALTY(1,2) 

F10.2     31-40    Mercury  (HG)  (yg/1)  QUALTY(1,3) 

F10.2     41-50    Suspended  sol  ids  (SS)  (mg/1)         QUALTY(1,4) 

F10.2     51-60    Volatile  suspended  solids  QUALTY(1,5) 

(VSS)  (mg/1) 

F10.2  61-70  BOD   (B0D5)  (mg/l)  QUALTY(1,6) 

F10.2  71-80  Total  organic  carbon  (T0C)(mg/l)      QUALTY(1,7) 

16  F10.2  1-10  COD  (COD)  (mg/1)  0JJALTY(1,8) 
F10.2  11-20  TKN  (TKN)  (mg/l)  QUALTY(1,9) 
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Table  4.   (continued) 


Card 

Card 

group 

Format 

col umns 

F10. 2 

21-30 

F10.2 

31-40 

F10.2 

41-50 

F10.2 

51-60 

F10.2 

61-70 

F10.2 

71-80 

17 

F10.2 

1-10 

F10.2 

11-20 

F10.2 

21-30 

F10.2 

31-40 

F10.2 

41-50 

F10.2 

51-60 

F10.2 

61-70 

F10.2 

71-80 

18 

F10.2 

1-10 

F10.2 

11-20 

F10.2 

71-80 

Descri  ption 


19 


20 


14 


A5 


1-4 


1-5 


N02  (N02)  (mg/1) 

NO   (N03)  (mg/1) 

Total  nitrogen  (TN)  (mg/1) 

TP0,  (TP04)  (mg/1) 

Chlorides  (CL-)  (mg/1) 

Lead  (PB)  (mg/1) 

Zinc  (ZN)  (mg/1) 

Tin  (SN)  (mg/1) 

Iron  (FE  (mg/1) 

Copper  (CU)  (mg/1) 

Cadmium  (CD)  (mg/1 ) 

Chromium  (CR)  (mg/1) 

Total  Sol  ids  (TS)  (mg/1) 

Oil  and  grease  (0&G)  (mg/1) 

PCB  (PCB)  (mg/1) 

Pesticides  (PEST)  (mg/1) 

Flow  (in  CFS)  at  the  time  this 
sample  was  taken  if  taken  at 
a  different  time  from  those 
points  input  in  card  group  M\ 

_99-,  this  value  terminates 
discrete  quality  input.  (One 
card  to  terminate  discrete 
quality  must  be  included.) 

Pollutograph  plot  option  in- 
dication ('PLUTO1  if  plot  is 
desired)   (if  I  OPT  =  'PLUTO' , 
include  card  groups  21  and  22) 


Variable 
name 

QUALTY 

[1,10) 

QUALTY 

[1,11) 

QUALTY 

[1,12) 

QUALTY I 

[1,13) 

QUALTY I 

[1,14) 

QUALTY I 

1,15) 

QUALTY 

[1,16) 

QUALTY I 

1,17) 

QUALTY 

[1,18) 

QUALTY < 

1,19) 

QUALTY 

[1,20) 

QUALTY I 

[1,21) 

QUALTY 

[1,22) 

QUALTY 

[1,23) 

QUALTY 

[1,24) 

QUALTY I 

'1,25) 

QUALTY I 

1,31) 

OPT 
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Table  4.   (continued) 


Card 

Card 

Group 

Format   Columns 

21 

1st  card 

20A^ 

2nd  card 

22 


23 


5A4 


1st 

card 

20A4 

2nd 

card 

5M 


Ik 


1-10 


Descr i  pt ion 


Discrete  quality  abbreviation 
from  card  group  15,  16,  17,  18 
(right-justified)  for  each 
pollutograph  to  be  plotted. 
Units  for  this  set  of  plots 
will  be  (mg/1).  Second  card 
must  be  included  even  if  blank 

Discrete  quality  abbreviation 
from  card  group  15,  16,  17,  18 
(right-justified)  for  each 
pollutograph  to  be  plotted. 
Units  for  this  set  of  plots 
will  be  (lb/min).  Second  card 
must  be  included  even  if  blank 

Sample  ID  no.  of  composite 
quality  (usually  an  ID  no. 
assigned  by  the  laboratory 
performing  the  quality 
analyses)  . 


Variable 
name 

INPUTQ 


INPUTQ 


SAMPNM 


E10.4 


F10.2 


11-20 


21-30 


F10.1  31-*K> 

F10.2  41-50 

F10.2  51-60 

F10.2  61-70 

F10.2  71-80 


Coliform  value  for  the  com-      QUALTY(l.l) 
posite  sample  (C0LF) 
(number/100  ml ) 

pH  value  for  the  composite       0_UALTY(1,2) 
sample  (pH) (A  value  must  be 
input  and  must  be  greater 
than  0.0) 

Mercury  (HG)  (ug/1)  QUALTY(1,3) 

Suspended  solids  (SS)  (mg/1)     QUALTY(1,4) 

Volatile  suspended  solids        QUALTY(1,5) 
(VSS)  (mg/1) 

BOD   (B0D5)  (mg/1)  QUALTY(1,6) 

Total  Organic  Carbon  (T0C)       QUALTY(1,7) 
(mg/1) 


ko 


Table  4.   (continued) 


Card 

Card 

Variable 

group 

Format 
F10.2 

col umns 
1-10 

Description 

name 

24 

COD  (mg/1) 

QUALTY(1,8) 

F10. 2 

11-20 

TKN  (mg/1) 

QUALTY(1,9) 

F10.2 

21-30 

N02  (N02)  (mg/1) 

QUALTY(1,10) 

F10.2 

31-40 

NO   (N03)  (mg/1) 

QUALTY(1,11) 

F10.2 

41-50 

Total  nitrogen  (TN)  (mg/1) 

QUALTY(1,12) 

F10.2 

51-60 

TPO^  (TP04)  (mg/1) 

QUALTY(1,13) 

F10.2 

61-70 

Chlorides  (CL-)  (mg/1) 

QUALTY(1,14) 

F10.2 

71-80 

Lead  (PB)  (mg/1) 

QUALTY(1,15) 

25 

F10.2 

1-10 

Zinc  (ZN)  (mg/1) 

QUALTY(1,16) 

F10.2 

11-20 

Tin  (SN)  (mg/1) 

0_UALTY(1,17) 

F10.2 

21-30 

Iron  (FE)  (mg/1) 

QUALTY(1,18) 

F10.2 

31-40 

Copper  (CU)  (mg/1 ) 

QUALTY(1,19) 

F10.2 

41-50 

Cadmi  urn  (CD)  (mg/1 ) 

QUALTY(1.20) 

F10.2 

51-60 

Chromium  (CR)  (mg/1) 

QUALTY(1,21) 

F10.2 

61-70 

Total  Solids  (TS)  (mg/1) 

QUALTY(1,22) 

F10.2 

71-80 

Oil  and  Grease  (OsG)  (mg/1) 

0_UALTY(1,23) 

26 

F10.2 

1-10 

PCB  (PCB)  (mg/1) 

QUALTY(1,24) 

F10.2 

11-20 

Pesticides  (PEST)  (mg/1) 

QUALTY(1,25) 

Include  blank  card  if  no  data 

27     F10.2      1-10    -99-  This  value  terminates 

the  composite  quality  input, 
and  the  job 
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Table  5-   Guidelines  for  executing  DSP. 


The  following  are  a  few  guidelines  for  executing  the  data  storage 
program. 

— Usually  only  rainfalls  greater  than  0.1  inches  (0.13  cm)  are   used 
when  calculating  the  number  of  dry  days  before  an  event. 

--The  start  of  precipitation  (card  group  3)  must  correspond  to  the  time 
of  rainfall  start  (card  group  11)  and  must  not  be  greater  than  the 
flow  start  time  (card  group  13). 

--Card  6  must  be  included;  it  should  be  blank. 

--The  flow  record  does  not  stop  with  the  end  of  the  precipitation  period, 
but  should  be  continued  for  a  reasonable  period  afterwards  to  ensure 
proper  total  volume  for  the  lb/min  composite  calculations. 

--Repeat  cards  15  through  18  for  each  time  interval ;a  discrete  sample 
was  taken  even  though  some  of  the  parameters  were  not  analyzed. 

--If  no  discrete  samples  were  analyzed,  only  composite,  delete  cards 
15  through  18  and  use  card  19  only. 

--Card  20  must  be  included  whether  pol lutographs  are  to  be  generated  or 
not. 

--If  the  pollutograph  plot  option  is  called,  cards  21  and  22  (both 
consisting  of  two  cards)  must  be  included. 

--If  an  oil  and  grease  grab  sample  is  taken  at  the  same  time  as  a 

discrete  sample  (two  separate  samples  with  different  sample  numbers), 
the  discrete  parameters  (four  cards)  must  precede  the  sample  number 
containing  only  the  oil  and  grease  (four  cards). 

--If  the  storm  event  under  consideration  is  a  snowmelt  (i.e.,  no  input 
hyetograph) ,  card  11  should  be  BLANK  and  card  12  should  have  a  -99- 
in  columns  1  to  10. 

--If  an  event  is  a  snowmelt,  do  not  use  the  water  equivalent  for  that 
melt  as  total  rainfall  (i.e.,  leave  columns  1-10  on  card  group  5 
BLANK) . 

— On  card  23  there  must  be  a  value  for  pH  and  it  must  be  greater  than 
0.0. 


hi 


An  example  of  the  input  data  for  an  individual  rainfall  runoff  event  is 
shown  in  Table  6.   Card  group  numbers  have  been  shown  in  the  right  hand 
column  as  a  means  of  identifying  input  data.   The  output  and  plots 
shown  later  in  this  section  will  have  this  data  as  input. 

DESCRIPTION  OF  THE  OUTPUT 


The  output  from  the  data  storage  program  is  largely  an  organized  review 
of  the  input  data.   Units  are  conveniently  shown  in  both  standard  and 
metric  systems.   Associated  sample  time  and  laboratory  number  are   dis- 
played alongside  the  data  for  easy  access.   The  output  is  listed  in  nine 
sections  namely: 


1 .  Title  block 

2.  Plotting  signal 

3.  Site  characteristics 
k.  Dustfal 1  analysis 

5.  Rainfal 1  record 


6.  Flow  record 

7.  Quality  parameters 

8.  Composite  parameters 

9.  End  of  run  signal 


Each  output  section  will  be  listed  on  at  least  one  individual  page, 
whether  there  is  input  data  to  be  displayed  or  not.   The  number  of  total 
output  pages  is  directly  proportionate  to  the  amount  of  input  data  and 
number  of  options  used. 

The  title  block  section,  Table  7,  identifies  the  program  source  as 
Department  of  Transportation  and  identifies  the  site  by  highway  name, 
snowmelt  event  (if  applicable),  date  of  storm,  beginning  of  precipita- 
tion (in  military  time)  and  a  unique  event  number.   All  but  the  first 
are    input  by  the  user.  The  program  source  name,  highway  name, 
snowmelt  event  (if  applicable)  and  event  number  are  also  displayed  at 
the  top  of  sections  two  through  eight. 

The  plotting  signal  section,  Table  8,  is  designed  to  be  displayed  if 
the  user  has  chosen  to  put  the  plotting  option  into  effect.   That  is, 
whether  the  user  intends  to  make  a  hyetograph,  hydrograph  and/or 
pol 1 utograph. 

The  site  characteristics,  Table  9,  ^^^    largely  user  input.   They 
include  storm  date,  precipitation  period,  start  of  overflow,  drainage 
area  in  acres,  total  precipitation  (for  this  storm  only),  duration 
of  precipitation,  number  of  dry  days  before  the  storm  and  number  of 
dry  days  to  an  accumulation  of  one  inch,  average  daily  traffic  during 
dry  days  and  average  daily  traffic  during  dry  days  to  an  accumulation 
of  less  than  one  inch  of  rainfall,  and  time  interval  between  discrete 
flow  samples  (in  minutes)  for  sampler  1  and  for  sampler  2.   The  program 
will  convert  and  display  total  precipitation  in  inches  and  centimeters. 
Total  rainfall  volume  over  the  entire  drainage  area    (in  cubic  feet)  is 
also  calculated  and  displayed.   It  should  be  noted  that  while  cubic 
feet  will  often  be  reported  in  six  or  more  digits,  it  is  only  significant 
to  the  number  of  digits  in  the  smaller  of  the  two  multipliers  (total 
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rainfall  and  site  area). 

The  dustfall  analysis  section,  Table  10,  prints  the  number  of  days 
the  dustfall  buckets  were  up  and  the  amount  of  sample  (in  mil il iters) 
and  dustfall  concentration  (in  miligrams  per  liter)  in  buckets  number 
1,  2  and  3-   All  of  the  above  are  user  supplied  input.   Total  grams 
are  calculated  for  bucket  number  1.   If  bucket  number  1  is  the  only 
bucket  operational  the  program  skips  to  the  end  of  the  dustfall 
analysis.   If  bucket  number  2  or  buckets  number  2  and  3  sre   operational 
the  program  will  repeat  this  above  sequence  and  then  go  to  the  end. 
At  the  end  of  the  dustfall  analysis  the  program  calculates  amounts  of 
dustfall  in  average  grams  and  grams  per  square  meter  per  day. 

The  rainfall  record  section,  Table  11,  displays  and  underlines  the 
interval  used  (in  minutes).   A  chart  with  four  columns,  time 
(military),  precipitation  (in  inches  and  centimeters),  intensity 
(in  inches  per  hour  and  centimeters  per  hour)  and  cumulative  precipita- 
tion (in  inches  and  centimeters)  is  then  displayed.   The  rainfall 
start  time  (military),  time  (interval)  and  precipitation  record  are 
user  supplied,  all  others  are  calculated  by  the  program. 

The  flow  record  section,  Table  12,  is  identified  by  the  heading  "flow" 
which  is  underlined.   Below  this,  three  columns  are  listed:   time 
(military),  cubic  feet  per  second  and  cubic  meters  per  minute.   The 
user  supplies  the  start  time,  end  time  and  interval  between  flow- 
points.   The  program  will  use  this  to  generate  the  flow  times  for  the 
list.   The  corresponding  cfs  values  are  recalled  from  the  input  and 
converted  to  cubic  meters  per  minute  and  both  columns  are  listed 
opposite  their  respective  times.   Finally,  the  cubic  feet  and  cubic 
meters  are  totalled  and  displayed. 

The  quality  parameter  section,  Table  13,  is  identified  by  the  heading 
"quality"  which  is  underlined.   Below  the  heading  are  several  columns. 
Time  (military)  and  sample  number  are   always  the  first  two  columns 
from  the  left.   Three  quality  parameters  can  be  displayed  in  columns 
to  the  right  of  this.   These  columns  include  the  parameter's  name 
and,  except  for  pH,  its  concentration  in  milligrams  per  liter  and 
pounds  per  minute.   If  more  than  three  parameters  are  input  a  new  row 
of  five  columns  will  be  displayed  below  the  first.   Again  the  first  two 
columns  from  the  left  are  time  (military)  and  sample  number.   Three  new 
parameters  can  then  fill  the  remaining  three  columns.   This  sequence  is 
repeated  until  all  input  quality  parameters  are  displayed.   If  a  value 
is  missing  (not  inputed)  a  ten  character  type  "#"  will  occupy  the 
otherwise  blank  space  in  the  column.   Time,  sample  number  and  milligrams 
per  liter  values  (or  pH  values)  for  each  quality  parameter  are  user 
suppl ied. 

The  composite  parameter  section,  Table  ]k,    is  identified  by  the  heading 
"composite"  which  is  underlined.   Below  are  up  to  thirty  quality 
parameters  in  five  rows  of  six  columns.   Each  displays  the  chemical 
analysis  information  about  its  respective  parameter.   The  parameter 
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abbreviations  appear  within  brackets  "<>".   User  supplied  concentrations 
are  shown  in  milligrams  per  liter,  micrograms  per  liter  for  mercury  or 
as  a  pH  value.   These  values  are  converted  (except  pH)  to  total 
kilograms  and  total  pounds  for  the  event.   The  sample  number  is 
printed  directly  below  the  heading. 

The  end  of  run  signal  section,  Table  15,  indicates  that  the  run  has 
completed  all  steps  of  the  data  storage  program.   That  is,  all  cards 
have  been  read  in  and  they  appear  to  conform  to  the  twenty-seven  card 
group  format. 

OPTIONS 

Highway  runoff  studies  by  their  very  nature  tend  to  generate  a  substan- 
tial amount  of  data  which  is  not  easily  handled  by  manual  procedures. 
The  Data  Storage  Program  has  been  set  up  to  utilize  a  Plotting  and/or 
Statistical  Option  to  capitalize  on  the  computational  potential  of  the 
already  stored  data. 

The  plotting  option  can  be  used  to  generate  reproduction  quality 
graphics  of  hyetographs,  hydrographs  and  pol lutographs  of  the  stored 
data.   This  option  is  particularly  useful  in  that  it  allows  an  analyst 
a  graphic  display  of  how  specific  monitored  parameters  varied  during  a 
rainfall-runoff  event. 

A  plot  storage  file  for  the  input  data  of  the  example  is  shown  in 
Table  16.   The  sequence  of  individual  plots  within  this  file  is  hyeto- 
graph,  hydrograph  and  pol 1 utographs  in  the  order  the  parameter 
abbreviations  were  assigned  on  cards  twenty-one  and  twenty-two.   The 
hyetograph,  Figure  7,  plots  rainfall  in  inches  per  hour  against  time 
(from  the  start  to  the  end  of  the  rainstorm).   The  hydrograph,  Figure  8, 
plots  flow  in  cubic  feet  per  second  against  time  (from  the  start  to  the 
end  of  the  flow).   The  pol 1 utographs ,  Figures  9  to  18,  plot  concentra- 
tions and  loadings  of  a  quality  parameter  in  milligrams  per  liter  and 
pounds  per  minute,  respectively,  against  discrete  sample  times  (from 
the  start  to  the  end  of  the  sampling  period).   There  is  a  possibility 
of  two  plots  per  each  parameter  when  plotting  pol 1 utographs--one  for 
concentrations  (mg/1)  and  one  for  loadings  (lbs/min).   Each  event  will 
consist  of  one  hyetograph,  one  hydrograph  and  selected  paired  pol luto- 
graphs.  That  is,  one  pair  for  each  parameter  set  inputted.   Any  one  or 
all  three  of  these  graphs  may  be  plotted  ijidsp-en-tly  of  the  others.   '>;  ', 
However,  once  the  storage  file  has  been  created,  the  pol 1 utographs  must 
be  plotted  in  the  sequence  defined  on  cards  twenty-one  and  twenty-two. 

In  our  example,  if  mercury  had  no  plot  file/,  our  first  plot  would  be 
suspended  solids  since  coliforms  and  pH  are  considered  meaningless  as 
plotted  data.   Plots  can  be  skipped  over  in  the  storage  file  until  a 
desirable  parameter  is  found.   However,  once  skipped,  the  file  cannot 
be  recalled  for  plotting  at  this  time.   The  event  would  have  to  be 
entirely  rerun  to  plot  a  skipped  pollutograph  plot  file.   Concentration 
plots  precede  loading  plots.   The  plots  illustrated  in  Figures  7  through 
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Table  16.   Example  of  plotting  file. 


1 

'LOTTING  DATA  FROM _DSP 

1 

MILWAUKEE  794 

24  617771910 

0 

13   13 

1910 

.0 

1915 

.1200 

1920 

.1200 

2 

5   16 

1925 

2.400 

5 

1930 

2.040 

1925 

.2100 

1935 

.8400 

1932 

.6400 

1940 

.8400 

_   1934 

6.600   .  ... 

1945 

.2400 

1939 

1.420 

1950 

.2400 

2004 

.2200 

1955 

.2400 

2000 

.1200 

2 

5   18 

2005 

.1200 

5 

2010 

.0 

1925 

.7700 

1932 

10.90 

1 

22  22 

1934 

82.40 

1925 

.1000E-01 

1939 

36.30  ...  .. 

1930 

.1000 

2004 

4.830 

1935 

3.650 

1940 

4.050 

3 

5   .  4 

1945 

2.400 

5 

1950 

1.650 

1925 

.4308E-03 

1955 

1.150 

1932 

.9800 

2000 

.7500 

1934 

19.83 

2005 

.5500 

1939 

7.250 

2010 

.5500 

2004 

.1939 

2015 

.3000 

2020 

.2000 

3 

5  22 

2025 

.1000 

5 

2030 

.1000E-01 

1925 

.1847E-01 

2035 

.1000E-01 

1932 

1.055 

2040 

.1000E-01 

1934 

45.97 

2045 

.1000E-01 

1939 

21.51 

2050 

.1000E-01 

2004 

.7080 

2055 

.1000E-01 

2100 

.1000E-01 

3 

5   15 

2105 

.1000E-01 

5 

2110 

.1000E-01 

1925 

.7492E-05 

1932 

.5244E-02 

2 

5    4 

1934 

.5462 

5 

1939 

.1521 

1925 

23.00 

2004 

.3933E-02 

1932 

654.0 

1934 

1307. 

3 

5   16 

1939 

391.0 

5 

2004 

69.00 

1925 

.3933E-05 

1932 

.9590E-03 

2 

5   22 

1934 

.1001 

5 

1939 

.2633E-01 

1925 

986.0 

2004 

.6181E-03 

1932 

704.0 

1934 

3030. 

3 

5   18 

1S»39 

1160. 

5 

2004 

252.0 

1925 

.1442E-04 

1932 

.1633E-01 

2 

5   15 

1934 

1.250 

5 

1939 

.6731 

1925 

.4000 

2004 

.1357E-01 

1932 

3.500 

1934 

36.00 

-999 

0 

1939 

8.200 

2004 

1  .  '<  0  0 
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18  were  made  using  a  TEXTRONIX  4014-1  catnode  ray  tube  and  printed 
with  a  TEXTRONIX  HARDCOPIER.   However  ."^mosTTg  r  a  phi  c  printers  can  be 
utilized  to  plot  data  generated  by  DSP. 

In  the  statistical  analysis  option  several  canned'  stati stical  packages 
are  available  including  the  Statistical  Analysis  System  (SAS) , 
Statistical  Package  for  the  Social  Sciences  (SPSS)  and  Biomedical 
Computer  Programs  (BMD) .   These  statistical  packages  contain  a 
wide  assortment  of  subroutines  which  will  perform  Analysis  of  Variance, 
Scatter  Plotting  and  Prediction  of  a  Best  Variable  Model,  to  name 
only  a  few  of  their  features.   An  example  of  the  statistical  analysis 
file  and  the  statistical  output  are  presented  in  Table  17  and  18,   /( 
respectively. 

The  data  storage  program  is  an  invaluable  tool  for  the  organization, 
conversion  and  correlation  of  raw  data  into  a  useful  and  accessible 
storage  unit,  from  which  unlimited  options  become  available  for  the 
prediction  of  quality  and  quantity  characteristics  associated  with 
highway  runoff. 
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Table  17-   Example  of  statistical  analysis  file, 


STATISTICAL  DATA*  REFER  TO  MANUAL  FOR  SPECIFIC  PARAMETERS 

MILWAUKEE  794.        24 
1.549 

1.430  __..  ...  ._ 

6.096 
.0 

.0  _ 

.0 

.0 

.0 

.3451 

132.1 

105.0  ..._ .•....;. ; 

1.201 

6.881 

15.00 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0  .  ...  ....  _ _..... 

6.800 

.S000E-03 

169.0 

144.0 

.0 

18.00 

94.00 

1.500 

.0 

.0 

.8400 

.4400 

30.00 

2.100 

.3600    __    , 

.0 

8.200 

.9000E-01 

.5000E-01 

.1000E-01 

407.0 

.0 

.0 
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FEDERALLY  COORDINATED  PROGRAM  (FCP)  OF  HIGHWAY 
RESEARCH  AND  DEVELOPMENT 


The  Offices  of  Research  and  Development  (R&D)  of 
the  Federal  Highway  Administration  (FHWA)  are 
responsible  for  a  broad  program  of  staff  and  contract 
research  and  development  and  a  Federal-aid 
program,  conducted  by  or  through  the  State  highway 
transportation  agencies,  that  includes  the  Highway 
Planning  and  Research  (HP&R)  program  and  the 
National  Cooperative  Highway  Research  Program 
(NCHRP)  managed  by  the  Transportation  Research 
Board.  The  FCP  is  a  carefully  selected  group  of  proj- 
ects that  uses  research  and  development  resources  to 
obtain  timely  solutions  to  urgent  national  highway 
engineering  problems.* 

The  diagonal  double  stripe  on  the  cover  of  this  report 
represents  a  highway  and  is  color-coded  to  identify 
the  FCP  category  that  the  report  falls  under.  A  red 
stripe  is  used  for  category  1,  dark  blue  for  category  2, 
light  blue  for  category  3,  brown  for  category  4,  gray 
for  category  5,  green  for  categories  6  and  7,  and  an 
orange  stripe  identifies  category  0. 

FCP  Category  Descriptions 

1.  Improved  Highway  Design  and  Operation 
for  Safety 

Safety  R&D  addresses  problems  associated  with 
the  responsibilities  of  the  FHWA  under  the 
Highway  Safety  Act  and  includes  investigation  of 
appropriate  design  standards,  roadside  hardware, 
signing,  and  physical  and  scientific  data  for  the 
formulation  of  improved  safety  regulations. 

2.  Reduction  of  Traffic  Congestion,  and 
Improved  Operational  Efficiency 

Traffic  R&D  is  concerned  with  increasing  the 
operational  efficiency  of  existing  highways  by 
advancing  technology,  by  improving  designs  for 
existing  as  well  as  new  facilities,  and  by  balancing 
the  demand-capacity  relationship  through  traffic 
management  techniques  such  as  buo  and  carpool 
preferential  treatment,  motorist  information,  and 
rerouting  of  traffic. 

3.  Environmental  Considerations  in  Highway 
Design,  Location,  Construction,  and  Opera- 
tion 

Environmental  R&D  is  directed  toward  identify- 
ing and  evaluating  highway  elements  that  affect 


*  The  complete  seven-volume  official  statement  of  the  FCP  is  available  from 
the  National  Technical  Information  Service,  Springfield,  Va.  22161.  Single 
copies  of  the  introductory  volume  are  available  without  charge  from  Program 
Analysis  (HRD-3),  Offices  of  Research  and  Development,  Federal  Highway 
Administration,  Washington,  D.C.  20590. 


the  quality  of  the  human  environment.  The  goals 
are  reduction  of  adverse  highway  and  traffic 
impacts,  and  protection  and  enhancement  of  the 
environment. 

4.  Improved  Materials  Utilization  and 
Durability 

Materials  R&D  is  concerned  with  expanding  the 
knowledge  and  technology  of  materials  properties, 
using  available  natural  materials,  improving  struc- 
tural foundation  materials,  recycling  highway 
materials,  converting  industrial  wastes  into  useful 
highway  products,  developing  extender  or 
substitute  materials  for  those  in  short  supply,  and 
developing  more  rapid  and  reliable  testing 
procedures.  The  goals  are  lower  highway  con- 
struction costs  and  extended  maintenance-free 
operation. 

5.  Improved  Design  to  Reduce  Costs,  Extend 
Life  Expectancy,  and  Insure  Structural 
Safety 

Structural  R&D  is  concerned  with  furthering  the 
latest  technological  advances  in  structural  and 
hydraulic  designs,  fabrication  processes,  and 
construction  techniques  to  provide  safe,  efficient 
highways  at  reasonable  costs. 

6.  Improved  Technology  for  Highway 
Construction 

This  category  is  concerned  with  the  research, 
development,  and  implementation  of  highway 
construction  technology  to  increase  productivity, 
reduce  energy  consumption,  conserve  dwindling 
resources,  and  reduce  costs  while  improving  the 
quality  and  methods  of  construction. 

7.  Improved  Technology  for  Highway 
Maintenance 

This  category  addresses  problems  in  preserving 
the  Nation's  highways  and  includes  activities  in 
physical  maintenance,,  traffic  services,  manage- 
ment, and  equipment.  The  goal  is  to  maximize 
operational  efficiency  and  safety  to  the  traveling 
public  while  conserving  resources. 

0.  Other  New  Studies 

This  category,  not  included  in  the  seven-volume 
official  statement  of  the  FCP,  is  concerned  with 
HP&R  and  NCHRP  studies  not  specifically  related 
to  FCP  projects.  These  studies  involve  R&D 
support  of  other  FHWA  program  office  research. 


DOT  LIBRARY 


00057CH3 


